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^ ■ Abstract 

. In the context of neutrino factories, we review the solution of the degeneracies in the neutrino 

oscillation parameters. In particular, we have set limits to sin'^ 2^13 in order to accomplish the 
I unambiguous determination of 823 and 5. We have performed two different analysis. In the first, 

psj ' at a baseline of 3000 km, we simulate a measurement of the channels Ve i^^, t'e — > i^r and 

, Un Un, combined with their respective conjugate ones, with a muon energy of 50 GeV and a 



running time of five years. In the second, we merge the simulated data obtained at L = 3000 km 
with the measurement of z^e — > i^fj. channel at 7250 km, the so called 'magic baseline'. In both 
cases, we have studied the impact of varying the i^r detector efficiency-mass product, {eu^ x Mr), 



QhI at 3000 km, keeping unchanged the detector mass and its efficiency. At L = 3000 km, we found 

the existance of degenerate zones, that corresponds to values of ^13, which are equal or almost 



equal to the true ones. These zones are extremely difficult to discard, even when we increase the 



^ , number of events. However, in the second scenario, this difficulty is overcomed, demostrating the 

Cd I relevance of the 'magic baseline'. From this scenario, the best limits of sin 2^13, reached at 3a, for 

sin^ 26*23 = 0.95 , 0.975 and 0.99 are: 0.008, 0.015 and 0.045, respectively, obtained at 5 = 0, and 
considering {e^^ x Mr) ~ 125, which is five times the initial efficiency-mass combination. 
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I. INTRODUCTION 



Nowadays, the neutrino oscillation is a confirmed phenomenon, supported by a huge 
amount of experimental data . More in detail, the neutrino oscillation frame- 

work is characterized by three mixing angles (^^12, ^23) ^13)) a CP-violation phase (6) and 
the squared mass differences (Am^i, Amg^). The mixing angles and the phase describe the 
PMNS mixing matrix, which is responsible for the connection between the mass eigenstates 
and the flavor eigenstates, while the squared mass differences play the role of frequencies in 
the oscillation |6|. 

Our present knowledge of all of these parameters goes as follows: the solar [l| and reac- 



tor neutrino experiments have obtained measurements of 612 and Am2i- Similarly, exper- 
iments with atmospheric and accelerator neutrinos j^, 4] have determined 623 and |Am|^|. 
Regarding ^13 or 6, no experiment has been able to measure any of these two parameters, 
even though an upper bound has been given to 6*13 Isl. 

Q n rl 

Common goals for present [T\ and future [8|, 19| neutrino experiments are the determi- 
nation of 6'i3, 6, the sign of Am|j^, as well as the precise measurement of the remaining 
parameters, in particular, whether sin^ 2^23 is maximal or not. However, it turns out that 
three types of ambiguities can appear when measuring oscillation parameters 10|: (^13, 5), 
sgnlAml^) and {623,71/2 — 623)- The combination of all ambiguities can produce an eight- 
fold degeneracy which impedes the unique determination of the involved parameters, and 
requires a combination of measurements to be realized in order to break the ambiguity. A 
big number of proposals has been taken into account since the degeneracy problem was 



identified, and many aspects o 



specific scenarios 
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In a long term view the experimental proposal with the best average sensitivity for 
measuring the oscillation parameters and unraveling their degeneracies is the neutrino fac- 
tory [29^ . These facilities, based on muon decay, will be able to produce a very high neutrino 
flux, and at the same time have the potential to measure six oscillation channels per polarity. 
Of these oscillation channels, the z/g 
measure 6*13 and S 



channel was identified as the one better suited to 



while z/^ gives very precise information of ^23 [iSj] . Besides, it is 
well known that a combination of these channels (e.g. u^. —>■ and v^. — > v^.) will be crucial 
for disentangling the degeneracies, especially helping us in the case of ^^23 is very close to 
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7r/4 (l0|. 

One of the most quoted distances for the location of the neutrino factory is 3000 km, 
due to its sensitivity when measuring S 30| and the sign of Am|]^ [3]]]. In fact, there are 



many references which discuss the degeneracy problem and its solution at this baseline 
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17|,ll8|,M. In this study, we review the solution of the degeneracy problem 
at 3000 km, using the neutrino factories and analyzing the channels z/g — >■ z/^, z/e ^ and 
Vfj_ —>■ as well as their conjugates channels. In our work, we point out the existence 
of degenerate zones that are very difficult to eradicate even when the number of events 
are increased. Moreover, it is shown how we can solve these problematic degenerate zones 
by the introduction of a second z/g Vu measurement at 7250 km, the so-called 'magic 



baseline' 



32] , and its combination with the simulated data obtained for 3000 km. 



A key feature in our analysis is the magnification of the number of events, in order to set 
limits on sin^ 2^13 to eliminate all degenerate solutions. This change in the number of events 
will be achieved through the gradual upgrading of the detector, this means that we will 
have a variable averaged detection efficiency and detector mass for v-r neutrinos, with initial 
values of ~ 5% and 5 Kton, correspondingly. For detection, the detector mass is fixed at 
50 Kton and the efficiency depends on the channel and is described in the Appendix. The 
muon energy is 50 GeV and the exposure time is five years. 

This paper is organized as follows: in section [TTl we present a brief reminder of the 
oscillation formalism and the probabilities of our interest. In section HITllIVI we present the 
solution of the degeneracies obtained from the analysis of the probabilities. Additionally, 
using only the probabilities, we recreate, as far as possible, the experimental conditions that 
we are going to use. This discussion will be very useful for understanding our experimental 
simulations. In section |Vl we will present our experimental simulation. Finally, in section |VT] 
and in the Appendix we present, respectively, our conclusions and all the experimental details 
considered in this work. 



II. THE NEUTRINO OSCILLATION PROBABILITY 

Neutrino oscillations are due to the non-coincidence between neutrino flavor eigenstates, 
Ua (a = e, /x,r), and neutrino mass eigenstates, z/j {i = 1,2,3, with mass mi, 1712, ms). 
These eigenstates are related by the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix U, 
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through: 



(1) 



C12C13 



■5l2Cl3 



"■S12C23 — Ci2S23Si3e*'^ C12C23 — Si2S23Si3e*'^ S23C13 



(2) 



where U is written using the standard parametrization 
/ 

U = 

^ S12S23 - 0120235136"" -C12S23 - ^120235136"" C23C13 

where Sij = sinOij and Qj = cosOij. We can see that U depends on three mixing angles, 6*12, 
6*23 and 6*13, and a CP phase, 6. Throughout this work, unless otherwise noted, we will use 
the following mixing parameters as 'true' parameters: 



(3) 



We shall work with the tangents of ^^23 and ^13, their 'true' values being: 

tan^ 023 = 6.3 x IQ-^ 
tan^ ^13 = 2.6 x IQ-^ 

In order to describe the evolution of neutrinos through matter, we use the evolution 
equation: 



sin^ 


26*12 — 


0.8 


sin^ 


26*23 = 


0.95 


sin^ 


2^13 = 


0.1 




6 = 


0°. 



(4) 



. d 

dx 



1 



( 



U 




Am^^ 









\ 



^ A{x) 0^ 



\ 








(5) 



where x = ct, Am?- 
by the potential 



and Ei^ is the neutrino energy. Matter effects are introduced 



A{x) = 2V2GfN^{x)E^ 

= 1.52 X 10~^N^{x)E^, 



(6) 



where Ne{x) is the electron number density. The 'true' values for Amfj we use throughout 
this work are: 

8.0 X 10"^ 
2.2 X 10"3, 



Am|i 



(7) 



following a normal neutrino mass hierarchy. 

Our theoretical analysis will use approximate formulae for neutrino oscillation probabil- 



ities derived in 



331] and expanded in 10|, l3J]. We use the notation of 10|. These formulae 



are valid in matter of constant density, where the average values of N^. are taken from lOj 



The formulae in 



id . I33I . [SJ] are expansions in 6*13 and a, where 



a 



If we use the following notation, 

A 

A = 
X = 

y = 
f = 

9 = 



= 1.27 



Am^^(eV^)L(km) 
KiGeV) 



A 



S23 sin 2013 

ac23 sin 26^12 
sm ( [1- A] A 



1-A 



sm I AAj 



A 



(8) 



(9) 
(10) 

(11) 
(12) 

(13) 
(14) 



where L is the distance between the neutrino source and detector, we can write the following 
oscillation probabilities: 



Pv^^u^ = (xf Y + {yg^ + 2xyfg (cos 5 cos A + sin 5 sin A) , 
Pue--ru^ = cot^ 02z{xfY + tan^ 02z{ygY - '^^yfg (cos 5 cos A + sin 5 sin A) 

P^^^^^ = sin^ 26*23 sin^ A + ^" ^ ^^3 sin^ 26*23 ^ / sin A cos AA - ^ A sin 2A 

+ a (2si3fg sin 2^12 sin 26^23 sin 6 sin A — c^2 sin^ 2^23^ sin 2A) 

— a— Si3 sin 2^12 sin 26*23 cos 26*23 cos 6 sin A (a sin A — g cos(yi — 1) A 

A-1 ^ 



(15a) 
(15b) 



.,2 ^^^2 I 

2A' 



— a" — - sin 2^12 sin^ 2^23 

g sin A cos(A — 1) A sin 2A 



A 
~2 



+ a^c\^ sin^ 2^23 cos 2 A 



(15c) 



The corresponding expressions for antineutrinos are equivalent to Eqs. (fTSll . replacing 
A —>■ —A and 6 —>■ —6. For a time reversed channel one needs to replace 6 —6. All 



formulae are valid for both normal and inverted neutrino mass hierarchies. Note that us- 
ing combinations of these three channels, one can calculate the probability of any other 
oscillation channel. 



III. DEGENERACIES IN NEUTRINO OSCILLATION PARAMETERS 

Given an oscillation probability function for the channel z/q, ^ z/^, defined by P^^^i^^^ = 
Fai3 {Oij, Amli, 6; L, E^), we have a degenerate solution if, for a given value P': 

P' = F„^(...,r/,C,...;L,E,) = F,^(...,r/',C',---;^,^.), (16) 

where (77, () and [i]', (') are different values for the same oscillation parameters. This means 
that there may be several different sets of oscillation parameters that could satisfy the same 
'measured' value of probability, being thus unable to distinguish the real set. 

It is important to note that a degeneracy is defined within the same oscillation channel, 
and for the same values of L and E,^ (experimental configuration). This indicates that the 
only way to solve any present degeneracy is by combining either several oscillation channels 
or different experimental configurations. 

There are three degeneracies associated with the standard — » i/^ channel measurement 



at neutrino factories: (6^13, 5), sgn^Amlj) and (6*23, 7r/2 — 6*23) [10 1. 



Having a (^13, 6) degeneracy, also called 'intrinsic' degeneracy, means having different sets 



of parameters ^13 and 6 giving the same oscillation probability [20|]. In particular, a CP- 
conserving scenario (sin 5 = 0°) can have a CP-violating (sin 5 7^ 0°) degenerate solution. 
As both parameters appear always together in the PMNS matrix, having a small ^13 can 
make the measurement of 6 much more difficult. 

The sgn{Am'^i) degeneracy rises due to our ignorance of the neutrino mass hierarchy, 
where we can have a normal {m^ > mi) or an inverted hierarchy {m^ < mi). This can lead 
us into two possible sets of (^13, 6), one for each sign of Am|i, which, when combined with 



the intrinsic degeneracy, results in a four-fold degeneracy 351] 



The {623,71/2 — 623) degeneracy appears in the measurement of the z/^ channel. 



which can determine sin^ 26'23 |36| . The z/g measurement involves a much more complex 
degeneracy in 6^23 which might include (6*23, 7r/2 — 6*23) sets. However, z/^ — experiments 
will be carried out before neutrino factories are developed, reducing this degeneracy in 623 



6 



into the discrete {O2Z, vr/2 — ^23)- 

Since the real value of 62^ is very close to 7r/4, the (6^23, 7r/2 — 6*23) degeneracy could be 
the most difficult one to solve. In some cases one would need a very good experimental 
resolution in order to separate ^^23 from its degenerate value. 

Considering the three different degeneracies, a single measurement could include as much 
as an eight-fold degeneracy. This is shown in Figure [H 



This problem has been studied through different perspectives 
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12j, and there are 



currently many alternative procedures that attempt to solve one or more of these degen- 



eracies. These procedures suggest the use of matter effects [13, 
alternative channels such as Vp — ^ v^, v„ v., a nd /or z/^ 



;he measurement of 



combination of different baselines [20|, 121|, [23, 
ues of (J^v) within one superbeam experiment 



23 



15|, 



16 



17 



19|, the 



24l |. and the comparison of several val- 



atmospheric neutrinos 26 ^ 



251]. Other solutions contemplate the use of 



27( 1 , the addition of Beta-Beam and Superbeam data 



and a 



careful selection of EyjL [10|] in order to solve the degeneracies. 



IV. GENERATION AND ANALYSIS OF EQUIPROBABILITY CURVES 

Our analysis will concentrate initially on the dependence between tan^^is and tan^^23, 
obtained from the Vf, — ^ z/^ channel. Given a 'measured' probability P^^-^y^ = Pi, we can 
approximate sin 6*13 ^ ^13 and cos ^13 ~ 1, so that we can solve Eq. f llSal) for ^^13: 

= . [\\^^-{v'gY^xry\b-^) - y'gcos{6 - A) ) , (17) 

2/tan6'23 ^ V cos^ 6'23 J 

where we have made y = cos ^23 in order to make evident the dependence on ^23- 

Using this equation, we plot in Fig. [2] a curve in the parameter space tan^ 6*23 vs tan^ 613, 
called 'equiprobability curve'. In order to obtain this curve, we have taken Pi as the respec- 
tive probability for sin^ 2^13 = 0.1 and 0.015, with the other parameters fixed at their 'true' 
values, considering L = 3000 km and = 30 GeV. As we can see, there is a continuous 
range of values of 6'23 and ^^13 that satisfy Eq. f|T7|) . meaning that, for a given 6, we have a 
continuous degeneracy in (6^23, 6*13). It is obvious from Fig.[2]that this degeneracy can not be 
solved by combining neutrino and antineutrino channels, since both curves are very similar. 
It is possible to obtain a similar relation between tan^ ^23 and tan^ ^13 for the Uf, 
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channel ^. We proceed in analogous manner, solving in this case Eq. (115bp for ^13, and fixing 
Pu.^u^ = -P2, obtaining: 



In Fig. [3] we show the combination of equiprobability curves obtained from Eqs. (fT7|) and 
f|T8|) . Notice that both equiprobability curves intersect at one point, which determines the 
real value of 6'i3 and 6*23. Thus, apparently, the combination of z/g — > z/^ and z/g — > z^t- channels 
can solve this continuous degeneracy completely. However, this is not really the case, since 
these results were obtained assuming as known the values of 6 and s(7n(Am|]^), which are 
parameters linked to ambiguities. For that reason, we will vary each one of them to see how 
it will affect our results. We start with 6, closely related 613, to study its infiuence on the 
location of the intersection point. 

The set of intersection points corresponding to the variation of 6 is shown in Fig. HI Only 
the values of tan^ 623 that fall within the 3a range in [s?] are evaluated. The horizontal lines 
show us that all intersections take place for the same value of tan^ 613, no matter if 623 and 
6 are at their 'true' values or not. This means that even though 623 and 6 are not known, 
we can determine 6*13 accurately and decouple it from 6 by crossing — >■ i'^ and z/g — z/^- 
channels. 

A combination of measurements can also determine the sign of Am|]^. Using the fact 
that the same value of tan^ ^13 should be obtained with neutrino and antineutrino channels, 
we compare both channels assuming the two possible signs of Am|]^. For the results to be 
consistent, the lines corresponding to neutrinos and antineutrinos should overlap, indicating 
the same value of tan^ ^13. If this is not the case, the choice of sign is considered inadequate. 

In Fig. [5] we plot the intersection points for neutrinos and antineutrinos, considering both 
signs of Am|]^. We can see that the lines corresponding to the correct sign of Am|^(+) lie 
one on top of the other, while the position of the lines corresponding to the wrong sign of 
Am|^(— ) depend on the use of neutrinos or antineutrinos, giving different values of tan^ ^13. 
This result allows us to distinguish solutions using the correct and wrong sign of Aml^. 

So far, we have determined analytically tan^ ^13 and the sign of Arn^^, having a continuous 
degeneracy between tan^ 6*23 and 6. Once the former two have been determined, we shall see 



^ This channel has been identified in 



15| as a very useful tool for solving the (6*13 , S) degeneracy 
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that this continuous degeneracy can be solved analytically by combining different L or 
configurations. In particular, we shall concentrate on variations of Ey. 

In Fig. Owe plot, fixing ^13 and for neutrinos and antineutrinos, the dependence of tan^ 
on 5, for L = 3000 km and E,^ = 20, 30 and 40 GeV. These curves are obtained from the 
intersection points of z/g — > and z/g — i/^. We also show equiprobability curves for the 
{V^ — s> V^) channels for the same distance and Ey = 30 GeV. The values of sin^ 2^13 
used are 0.1, 0.05, 0.025 and 0.015. In order to differ the curves got from the intersection 
point of z/e — > z/^ and z/g from those of ^ v^-, we shall denote the former as v^®Vt-. 

The z/^ ® Vr equiprobability curve shows a continuous degeneracy, with any two neutrino 
(or antineutrino) curves at different energies intersecting at two points, a true and a fake, 
correspondingly. This means that each set of curves, corresponding to a particular ^13, has 
its own (6^23, 5) ambiguity, different from (6^23, 7r/2 — ^^23)- We will denote this degeneracy as 
(^23,'^)!^ and (6'23,5)i7 for neutrino and antineutrino channels, respectively. It is important 
to note that the bottom plots of Fig. [6] do not show a second intersection because it occurs 
outside the currently allowed values of tan^ 6*23 3^] ■ However, the fake points for neutrinos 
and antineutrinos are not the same, and they do not necessarily correspond to the degenerate 
value 71/2 — 62^1 marked by the z/^ z/^ channel. 

From the exposed above, we can infer that the minimal configuration for solving the de- 
generacies, analytically, consists in combining two different energies for v^^v^-, for neutrinos 
and antineutrinos. 

In Fig. El we show more extended information than this minimal configuration. This 
information will be useful for understanding our results of the simulation of the experimental 
setup. This is the case of the z/^ ® curves presented for three different energies for 
neutrinos and antineutrinos, which serve to mimic, roughly, a spectral analysis. Additionally, 
the introduction of z/^ channel is very important considering its statistical power to 

determine, precisely, 6'23 (even though, this is also valid for its degenerate value). The 
next section will have the purpose to discuss all of this information, visualized from an 
experimental perspective. 



9 



A. Experimental Considerations 

The previous section explains how to solve the degeneracies analytically. However, in an 
experimental context where statistical and systematic errors are present, the situation will 
be different. For example, whenever there are degenerate intersections of equiprobability 
curves, a low region around these points can appear during the fitting of data. Thus, it is 
possible to make rough predictions on how well the degeneracies will be solved by observing 
the confluence of equiprobability curves and their true and degenerate intersections. 

The kind of situation described above arises in Fig. [6l where the neutrino (antineutrino) 
z/^ z^T curves have a (6*23,5)1/ ((6'23j<^)f) degeneracy. Even though analytically these two 
degeneracies are not a problem, we can expect two low degenerate regions close to each 
other, one for each degeneracy. Fortunately, the inclusion into our analysis of the 
channel, with its large statistical weight, introduce the (^23? 7r/2— ^23) degeneracy, eliminating 
all of these (6*23, 5) degeneracies that do not coincide with (6*23, vr/2 — 6^23)- However, if any 
of the two ambiguities does coincide, one can expect a very low region which should be 
very hard to resolve. 

Another important issue, for having a perspective of the experimental behavior using the 
equiprobability plots, is the fact that not a unique value of 6*13 is going to be compatible 
with the experimental data. Therefore, it is necessary to generate ® curves for a range 
of values of 6*13. If this variation makes any (6*23,5) degeneracy overlap the (6*23, vr/2 — 6*23) 
degeneracy, ambiguous regions can appear in the fitting of data. From this variation of 6*13, 
we can expect a four-fold ambiguity to be present, which correspond to a pair of (6*13,5) 
degeneracies per 6*23 and 7r/2 — 6^23 ^• 

In Figures [71 to [TOl we show degeneracies for sin^ 26^13 = 0.1, 0.05, 0.025 and 0.015, re- 
spectively, the same values we used in Figure [61 These Figures contain four panels, each 
panel representing one of the four-fold ambiguities. We have generated the equiprobability 
curves using the probability corresponding to each given value of sin^2^i3. Keeping this 
probability as a constant, equiprobability curves for different (fake) values of sin^ 20^3 have 
been generated. The variation of this parameter shows four overlaps between intersection 
points, which represent the four-fold ambiguity. 

^ Notice that we can obtain an eight-fold ambiguity if we also perform this variation for the wrong sign of 
Am|]^. We will not do this, since this degeneracy will not be relevant in the considered cases. 
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In all cases, the upper left panel of each Figure shows the equiprobability curves using 
the real value of sin^ 26*13 (the true solution). Obviously, in this case, the z/^ z/^ curves 
overlap the common intersection of neutrino and antineutrino z/^ ® i^r curves, indicating the 
real values of tan^ 623 and 6. 

The upper right panel shows the displacement of the common intersection into 7r/2 — ^23, 
done by changing the value of sin^ 2^13 with which we generate the equiprobability curves. 
In this case, all ® Ur equiprobability curves intersect with each other on the degenerate 
6*23, with no noticeable variation in S, making the experimental resolution of this case quite 
difficult. This type of ambiguity shall be named 'type I mixed degeneracy', and is equivalent 

to {923,71/2-923). 

The lower panels of each Figure show the other two degeneracies. In these cases there is 
no common intersection between the neutrino and antineutrino z/^ (S> z/^ curves, but only the 
{923, degeneracy overlaps one of the z/^ z/^ curves or, in other words, either coincides 
with ^23 or 7r/2 — 6^23- Of course, this will be another case difficult to solve experimentally, 
worsened by the proximity of the (6*23, 5)i7 degeneracy. A similar situation will occur when 
(^'23, 5)f overlaps a.u^ ^ curve, however, since the value of sin^ 26*13 is close to the previous 
one, these results are not presented. 

In more detail, the lower left panel shows that the (6*23, 6)^ intersections lie on the — >■ 
equiprobability curve where 6*23 is at its real value. In this case, the degenerate 5 is noticeably 
different from the real one. We name these cases 'pure degeneracies'. 

The lower right panel shows a second value of sin^ 26*13 where the (6*23,5)1/ intersection 
overlaps the z/^ — > z/^ equiprobability curves, this time for 7r/2 — ^23- The value of 5 changes 
noticeably, like the previous case. These situations are named 'type II mixed degeneracies'. 

It is worthwhile to remark that Figs. [TlfTOl are going to be very useful in understanding the 
full experimental behavior that will be presented in Section |Vl In fact, the true values chosen 
to generate results will be equal to those we are analyzing in Figs. [TlfTOl In particular, we will 
be able to correlate these Figures directly with the tan^ 6^23 vs. 5 panel from its experimental 
counterparts (Figure [12] to [Hj) . 

In general, we must say that a key feature for disentangling the degeneracies is how large 
the difference between the true value of ^13 and the fake one is. Actually, in an experimental 
situation, we will need a high enough number of z/g z/^ and z/g — > v-r events, which can 
make this difference significant from the statistical point of view. Putting this in our context. 
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in Figure [TJ the degenerate sin^ 29[^ are 0.145, 0.0668 and 0.09 for the 'pure,' 'type I mixed' 
and 'type II mixed' degeneracies. These correspond to a difference of about 45%, 33% and 
10%, respectively. Following our reasoning at the begining of the paragraph, within an 
experimental framework we can expect the 'pure' degeneracy to be the easiest one to solve 
(45% variation), while the 'type II mixed' should be the hardest one (10%). The order in 
the differences between the true value of 6*13 and the fake one, for each type of degeneracy, 
can be explained qualitatively, seeing, for example. Figure [3 There, we can observe that 
the 'type II mixed' (lower right) degeneracy is exhibiting the more alike behaviour in the 
pattern of the curves and position of the intersections showed by the true solution (upper 
left). In the case of 'type I mixed' degeneracy, which is the next in order, 6 coincides with 
the real value, an all the equiprobability curves intersect at n/2 — 623. The last in this list 
is the 'pure' degeneracy, where 6 is rather different respective to the real one, and only one 
type of equiprobability curve (neutrino or antineutrino) intersects at 6^23. 

An interesting case to keep on mind is shown in Figure [HI for sin^2^^i3 = 0.05. In this 
case the 'type II mixed' degeneracy occurs when sin^ 26*13 acquires its real value, the upper 
left and lower right panels are identical. In fact, we see that the true solution and the 'type 
II mixed' degeneracy are appearing in these panels. This means that our results will be 
ambiguous even if ^13 is known exactly. We call this situation a 'persistent' degeneracy, and 
happens only for certain true values of sin^ 2^13. If the real 6*13 is such that a 'persistent' 
degeneracy is present, we can expect the experimental fitting of data to show an ambiguity 
even when the number of events in all channels is large. To disentangle this problem, we 
will need an extra measurement made at a different baseline. 



1. The Magic Baseline 



Previous studies in [10|, show that there are certain baselines where the approximate 
expressions for the oscillation probability in matter f|T5l) can be greatly simplified. By 



choosing L such that g is zero, one finds for the z/g 



channel: 



(19) 



It was noted in [32'] that the 
liminary Reference Earth Model 



appropriate L for an experiment using the realistic Pre- 
is L = 7250 km. By doing this, 6*13 can be very well 
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determined, even though the sensitivity to S is lost. Thus, in order to measure all param- 
eters adequately, this experiment carried out at L = 7250 km must be complemented with 
one at L = 3000 km. The former has received the name 'magic baseline', due to its good 
capabilities in the measurement of and ^^23- 

Figure [TT] shows equiprobability curves for an experimental configuration corresponding to 
the 'magic baseline', in contrast to those for L = 3000 km, which are drawn in black. Curves 
are presented for sin^2^i3 = 0.1, 0.05, 0.025 and 0.015, but, as the analytic formulae f|T5|) 
are only valid for L < 4000 km [lo|, we generate these curves numerically. 




The curves from the 'magic baseline' clearly can not determine the value of 6, but are 
very useful in constraining the value of tan^ 623 within the proper octant. This is a very 
attractive characteristic when trying to solve the 'persistent' degeneracies, since it is only 
present for the degenerate 11/2 — 623. However, the number of events coming from this choice 
of baseline is not expected to be as significant in helping us to obtain a precise measurement 
of 623, as the number of events in the z/^ — z/^ channel at L = 3000 km, so it should not 
become a substitute for the latter. 

V. EXPERIMENTAL TESTS OF THE DEGENERACIES AT NEUTRINO FAC- 
TORIES 

In this section we are going to simulate the combination of oscillation channels and their 
corresponding baselines we proposed in preceding section, within the experimental context 
of a Neutrino Factory. 

The determination of the oscillation parameters is done through the measurement of 
an experimental observable: the number of events. This observable is a reflection of the 
transition probability convoluted with the different experimental factors, such as neutrino 
flux and cross-section. In this work we shall concentrate on solving the degeneracies within a 
neutrino factory context, following the procedure described in the previous section. Details 
of the generation of number of events can be found in the Appendix. 
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A. Description of Statistical Analysis 



The following shall be described in terms of neutrino fluxes produced from the decay of 
/i"*". Nonetheless, the analysis is also performed using neutrino fluxes from fi~ decay, and 
these will be implicitly included. 

We are going to present two types of analysis, one for the L = 3000 km baseline, and 
another one for a combination of L = 3000 km and L = 7250 km baselines. At the L = 
3000 km baseline, we assume the installation has two detectors, a magnetic iron calorimeter 
for the Up — >■ Uj, a nd — > z7^ channels and an emulsion cloud chamber for the Ue Vt 



channel, as in 39|]. At the L = 7250 km baseline, we consider an identical magnetic iron 



calorimeter measuring only z/g oscillations. 

To analyze the generated data at L = 3000 km, we follow the procedure outlined in [l^, 
HQ. For the i'^ and z/^ channels we use 20 energy bins with a width of 



2.3 GeV each. For the ^ channel we use four bins of variable width, as in [171]. For the 
L = 7250 km baseline, we do not perform any binning, and use the total number of events. 

Denoting and as the number of signal and background events generated for the 
channel — > 1^/3 in the ith bin, we perform either a Gaussian or Poissonian smear on each 
bin, depending on the number of events, to simulate statistical errors: 



)obs = Smear (n^^ + 6^^) (20) 



where (^^^)o6s is the 'measured' number of events for the channel z/^ in the ith bin. 

This 'measured' value is fit to the theoretical (non-smeared) number of events {n\fj)thi as a 
function of tan^ ^13, tan^ ^23 and 5. We either use a Gaussian or Poissonian minimization, 
depending on the value iji'^pjobs'- 



(21) 



1/ ("1/3)068 

^ 2 ({Kp)th - {Kp)obs + {Kp)obs In ^f^) {Kp)obs < 5 
In the first analysis the total is defined by: 

= E + + ix%,^v, + (^ - V)] (22) 
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and for the combination the is: 




J L=3000 km 



i 



+ kx')t%, + - 

L ° J L=7250 km 



(23) 



where we use the total number of events for the 7250 km basehne. 
B. Results for L = 3000 km 

Initial results can be seen in Figure [12], which were generated using an average global 
detection efficiency of eo ~ 5% and a detector mass of Mq = 5 Kton (see Ap- 
pendix IA3I) . Fig. [12] consists of four parts corresponding, from left to right and top to 
bottom, to sin^ 26*13 = 0.1, 0.05, 0.025 and 0.015, respectively. Each part is composed by 
three panels, representing the three projections in the parameter space considered in our 
analysis: tan^ 6^13 vs. 5, tan^ 6*13 vs. tan^ 6^23 and tan^ 6*23 vs. b. In these panels we display 
the allowed regions for I0 (3o") confidence levels in red (blue). 

Within the evaluated sin^ 26*13, we do not get rid of the degenerate regions at a 3(j 
confidence level. In general, we observe that as long as we diminish sin^ 2^13, up to two 
degenerate regions can appear. In addition, with the decreasing of sin^ 2^13, each degenerate 
regions worsens, which means that they can include allowed regions with a confidence level 
under 2cr. 

We note, in the upper left part of Fig. [121 the appearance of a first region at 3(j, for 
sin^ 2^13 = 0.1. Recall from Section [IV Al that this degeneracy corresponds to a 'type II 
mixed' degeneracy. In our analysis of equiprobability curves, we expected this degeneracy 
to be the hardest one to solve, since the required variation of sin^ 2^13 was of 10%. If we 
decrease the 'true' value of sin^ 26*13 down to 0.05 (upper right), this degeneracy appears at 



A second degenerate region shows up when sin^ 2^13 = 0.025 (lower left), at 3o". Unlike the 
first degeneracy, this one is a 'type I mixed' degeneracy, and reaches 2a when sin^ 26*13 = 0.15 
(lower right). Notice that there are no 'pure' degeneracies present. 

The appearance of 'type I mixed' and 'type II mixed' degeneracies and their growing 
statistical compatibility with the real solution are due to the decreasing of the real value of 
sin^ 26*13. This implies a fewer number of v^. events which, jointly with the similar value of 



2a. 
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the real and degenerate sin^25i3, complicate our ability in statistically distinguishing these 
degenerate regions from the real one. In contrast, our results do not present any 'pure' 
degeneracies because the associated degenerate value of sin^ 26i^ has larger deviations with 
respect to the real one, making it easy to exclude at the evaluated confidence levels. 

The sgn{/\m\^) degeneracy does not present a problem in the evaluated cases. In Table |T] 
we show the minimum obtained using the correct and wrong signs of Am|^. In all cases 
the lowest assuming the wrong sign is much larger than any of the evaluated confidence 
levels defined by the overall lowest x^. Thus, no extra regions appear when analyzing the 
observed events using the wrong sign of Am|^ . This makes it possible to discard four of the 
eight possible solutions, reducing the problem into a four-fold degeneracy. 

Since the ambiguities are due to our inability to determine sin^ 26*13 accurately, we need 
a higher number of Vr events in order to solve the degeneracy problem. This can be 
achieved by increasing the detector efficiency e,^^ and/or the detector mass Mr. We are 
going to present our following results in terms of an efficiency-mass combination ratio, 
feM = {^vr ^ ^t) I (eo X Md), concentrating in particular in r^M = 3 and r^u = 5. We 
have to stress out that the functional form of the linear efficiency described in Appendix IA3I 
is preserved, as well as the background levels. 

In Figure [T3] we display the results using r^M = 3. For sin^ 26*13 = 0.1 (upper left) the 'type 
II mixed' degeneracy is solved. In contrast, for sin^ 26^13 = 0.05 (upper right) the situation 
does not improve respect to Figure [121 keeping the same confidence levels for the degenerate 
region. For sin^ 26*13 = 0.025 (lower left), the 'type I mixed' degeneracy disappears, reducing 
the 'type II mixed' degeneracy to 3a. In the case of sin^ 26*13 = 0.015(lower right), the 'type 
II mixed' degeneracy is solved, while the 'type I mixed' degeneracy was reduced to 3cr. Thus, 
the only degeneracy left at 2a occurs when sin^ 2^13 = 0.05. 

For r^M = 5, shown in Figure [TH the degeneracies are solved in all cases except for 
sin^ 2^13 = 0.05, which remains at 2a. The existence of this degenerate region is independent 
of the number of measured Ur events. This particular type of behavior occurs when the value 
of sin^ 26*13 for the degenerate solution is very similar (or equal) to the true value of this 
parameter. We have called this case a 'persistent' degeneracy, and has been discussed in 
Section IVAl 

In Figure [15] we show the values of sin^ 26*13 that present some type of degeneracy, as a 
function of r^M- The red (blue) regions correspond to ambiguities not solved at 2a {3a), 
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while white regions do not present any degeneracies up to 3a. This analysis is done for 
tan^^zs = 0.63, 0.73 and 0.82 (sin^ 2^23 = 0.95, 0.975 and 0.99), in the top, middle and 
bottom row, respectively, and for 5 = 0°, 60° and 90°, in the left, center and right columns, 
respectively. 

It is important to mention that we classify a case as non-degenerate if there is no more 
than one allowed region for 6'i3 and 6 within their respective parameter space, provided 623 
is confined within one octant. At the same time, we classify a case as degenerate in 623 if 
both octants are occupied, even if this is done by only one single region. Notice that this 
involves excluding maximal ^23- 

From Fig.[T5]we observe that, for a fixed ^23, the upper limits of the 2a and 3a degenerate 
regions decrease as 6 increases. For sin^ 26*23 = 0.95 and r^M = 1? the 3a upper limits for 
sin^ 26*13 are 0.11, 0.095 and 0.055 for 6 = 0°, 60° and 90°, respectively, while for r^M = 5 
they are 0.08, 0.06 and 0.025. For sin^ 26*23 = 0.975 and r^M = 1, the same 3a limits are 
0.19, 0.175 and 0.09, and for r^M = 5, these are 0.16, 0.125 and 0.055. The results for 
sin^ 26*23 = 0.99 do not present any upper limits. 

In contrast, for a fixed 6, the upper limits increase with sin^ 2^23, which is understood since 
values of sin^ 26*23 closer to unity are harder to disentangle from their respective degenerate 
pairs. This is clearly seen for sin^ 26*23 = 0.95, 0.975, being over the allowed limits in sin^ 26*13 
when sin^ 2^23 = 0.99. 

The plots of Figure [T5l corresponding to 5 = 0°, are also exhibiting a particular behavior. 
For a high enough value of r^M, we can have regions free of degeneracies at 3a (white 
regions) limited from above and below with regions that contain degeneracies under 3a and 
2a (blue and red regions) . The existence of the lower degenerate regions is straightforward to 
understand, because a low value of sin^ 26*13 implies a low number of Ur events, which worsens 
our ability to resolve ambiguities. This last fact is apparently contradicted by the presence of 
the blue and red regions above the white one, for higher values of sin^ 26*13. This contradiction 
appears because these regions belong to degenerate cases where the fake sin^ 26*13 is very 
similar to the corresponding real one (or even equal, this is what we called 'persistent' 
degeneracy). Therefore, our ability to statistically distinguish these degenerate situations 
from the real ones gets worse. These regions become narrower with higher values of r^M, given 
that an increase in the number of u^- events eliminate, progressively, the degeneracies at the 
required confidence levels. Examples of this narrowing happen when sin^ 2^13 = 0.025 and 
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0.015, shown in Figures [T2] to UM Similar arguments can be used to explain the intermediate 
blue zone for 6 = 60°. For 6 = 90°, no intermediate zone exists. 

Even though it is desireable to achieve a high value of r^^, especially if the true value 
of sin^ 29i3 liss within a white intermediate zone, this does not solve the red and blue 
degenerate zones above them. Considering that these zones are related to a 'persistent' 
degeneracy, where the fake sin^ 2^13 is very close or equal to the real one, the only way to 
solve this situation is by including additional information that precisely specifies the octant 
where 623 belongs. In the previous section we showed analytically that the precise definition 
of tan^ 6*23 can be obtained by using a second Ue measurement at the 'magic baseline'. 

C. Results for L = 3000 km L = 7250 km 

The presentation of these results starts at Figure [IHl which shows the same information 
as Figure [121 but includes a measurement realized at the 'magic baseline' into the analysis. 
It is evident that the inclusion of this baseline proves to be effective not only in solving the 
'persistent' degeneracy, but also in improving the unambiguous determination of the real 
solution. Note that even though t^m = 1, the top left, top right, and bottom left parts are 
completely degeneracy- free, leaving only the bottom right part (sin^ 26'i3 = 0.015) with a 
'type I mixed' degeneracy due to the small number of Ur events. 

Figure [17] is equivalent to Figure [15], with the additional information from the 'magic 
baseline'. Note that r^M still refers to improvements of the u^- detector at L = 3000 km. 
In all cases, the ambiguities related to 'persistent' degeneracies have been eliminated, and 
the lower degenerate zones due to the small number of Ur events have had their upper 
boundaries reduced. The behavior of the boundaries has become simple, decreasing with an 
increase of r^M- The impact of this improvement is easily appreciated in the cases where 
sin^26'23 = 0.99, where upper boundaries for the degenerate regions can be established. The 
most noticeable case occurs when sin^ 26'23 = 0.99 and 6 = 90°, where the once unsolvable 
case acquires the same simple behavior. 

The variation of 6 does not produce important changes in the behavior of the limits. 
Nonetheless, we can observe that the upper limits decrease as 5 ^ 0°, in opposition to the 
case without the 'magic baseline'. This change of behavior is due to the disappearance of 
the 'persistent' degeneracy, which in the former case brought extra degenerate regions for 
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low values of S. For sm^2^23, in contrast, the upper limits decrease depending on how far 
away this parameter is from unity, which is consistent with our previous analysis. Notice 
that a change of scale in sin^ 2^23 = 0.95 was necessary in order to show the results properly. 

Table [IT] shows the lowest values of sin^ 29i3 that solve the degeneracies using this con- 
figuration, considering the minimum r^M = 1 and maximum r^M = 5 efficiency-mass com- 
bination used in this work. For r^M = 1 at 2cr {3a), the best lowest value is 0.01, for 
sin2 2023 = 0.95, 5 = 60° (0.023, for sin2 2e23 = 0.95, 5 = 0°), while the worst is 0.06, for 
both sin2 2^23 = 0.99 and 5 = 0°,90° (0.105, for sin2 2^23 = 0.99, S = 90°). 

For r,M = 5 at 2(T (Sd), the best lowest value is 0.005, for sin2 2^23 = 0.95, 5 = 0° (0.008, 
for sin2 2^23 = 0.95, 5 = 0°), while the worst is 0.04, for sin^ 2^23 = 0.99 and 6 = 90° (0.07, 
for sin^ 2^23 = 0.99, S = 60°, 90°). 

In general, it is easier to solve the degeneracies when 5^0° and sin^ 26'23 is less com- 
patible with the unit. The eventual deviations from this tendency are due to the random 
smearing effect. 

VI. SUMMARY AND CONCLUSIONS 

In this work we reviewed the degeneracy problem, and studied its solution within the 
framework of a neutrino factory, extracting limits at 2cr and 3a in sin^ 2^13, for different 
values of 6*23 and 6. This was done using two experimental configurations with a muon 
energy of 50 GeV and five years of exposure time. The first one considered a baseline 
at L=3000 km involving the Uf, —>■ z/^, Uf, —>■ and channels and their conjugate 

probabilities, while the second one was a combination of the information obtained at L=3000 
km mixed with the z/g channel at the 'magic baseline', L=7250 km. For z/^ detection, 

we kept along the analysis a constant detector mass and efficiency. In contrast, we varied 
the product of the Uj. detector mass and the efficiency, {e,^^ x M^). 

In this paper, our experimental simulations were preceded by an analysis which used 
only equiprobability curves and took into account an uncertainty in the determination of ^13 
and the effect of the energy spectra. Within this framework, and considering L =3000 km 
with the channels involved in this analysis, we identified the four-fold degeneracy related to 
combinations of (6^13,5) and {623,71/2 — 623), naming them 'pure,' 'type I mixed' and 'type 
II mixed' degeneracies, where 'mixed' made reference to the ambiguous value of ^23 in the 
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(^23, 7r/2 — ^23) degeneracy. We also predicted the order of difficulty in solving experimentally 
these degeneracies, which was, from lesser to greater, 'pure,' 'type I mixed' and 'type II 
mixed'. 

Thus, a good determination of ^^13 is an essential objective when solving the ambiguities. 
The degeneracies with a fake ^13 similar to the real one would be harder to solve, since they 
are more compatible with the (8) i^r data. However, at L = 3000 km we found a particular 
'type II mixed' degeneracy that can arise for the true ^13, which we called a 'persistent' 
degeneracy. Since this ambiguity does not need the variation of ^13 to manifest itself, an 
increase in the Vr statistics will be insufficient to solve the problem. We showed that a second 
measurement of an oscillation channel was needed, and identified a z/g — z/^ measurement 
at the 'magic baseline' as a suitable choice for a second experiment. 

These results were confirmed by our experimental simulation. For the L=3000 km anal- 
ysis, we observed degenerate zones for intermediate values of sin^ 2^13, which does not dis- 
appear even for r^u — 5. This is very interesting, since a common thought is to relate the 
solution of degeneracies with the size of sin^ 2^13. Nevertheless, when we added the second 
measurement realized at the 'magic baseline', we managed not only to get rid of the 'per- 
sistent' degeneracy, as wc predicted before, but also improved the state of other degenerate 
situations. In fact, for t^m — 5, we got limits of sin^ 26*13, at 3cr, for sin^ 2^23 = 0.95, 0.975 
and 0.99 which were: 0.008, 0.015 and 0.045, respectively, obtained a,t 5 — 0. 
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APPENDIX A: SIMULATION OF NUMBER OF EVENTS 

In an experimental situation it is not the oscillation probability which is measured, but 
a certain type of event that occurs due to this probability. The experiment counts how 
many times this event occurs and then the statistical analysis compares this number with 
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theoretical predictions within parameter space. The simulation of the number of events shall 
take into account the neutrino flux, transition probability and cross-section, as well as the 
efficiency and resolution of the detector. 

Within this work, the measurement of number of events is done in the context of a 
neutrino factory. The neutrino beam in this experiment is produced by the decay of muons: 

^ ^ (Al) 

Being this a leptonic decay, the neutrino flux can be determined precisely, which provides 
very good statistics. One of the greatest advantages neutrino factories have in comparison 
with other experiments is the lack of contamination in the decay. Since the muon has only 
one decay channel, the amount of neutrinos of a given flavor expected in the absence of 
oscillations can be very well described. Another advantage is that this experiment produces 
large, similar fluxes of and V^. (or and Vf.)-, currently unavailable in other neutrino 
experiments. 

Studies of neutrino factory fluxes and detectors have taken place in several papers, such 



as 



17,^ 



4011 to calculate the number of events in a 



dn'r. 



dE' 



401 1 . We use the method described in 
neutrino factory. 

The differential event rate for the channel z/^ — > z//3, for a given interaction type (where 
<^ = CC, NC), is given by: 

l = N JdE.x (l)^{E,) X P,^^,,{E,) X a^iE,) x Rf,{E,, E',) x e^^El) (A2) 

where E^, is the incident neutrino energy, E^ is the detected neutrino energy, is a nor- 
malization constant, (p^ is the flux at the detector, Pj^^^j^^ is the neutrino oscillation 
probability between a and (3 flavors, o"^ is the cross-section for the interaction type 
and i?/3 and are the energy resolution and detection efficiency, respectively, for the z/^. 

As a signal we will consider the z/g ^ t'^, z/g ^ z/^ and z7^ channels, and the 

respective channels for antiparticles, measured through charged current interactions. We 
also calculate the appropriate backgrounds using Eq. flA2p . replacing the efficiency with 
a background level. In the next sections we shall describe each of the factors that appear in 
Eq. (IA21). 
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1. Normalization Constant 



The normalization constant N includes experimental information that can be factorized 
out of the integral. In the case of neutrino factories, we use: 

N = Uy X 10^ X Na y< Md, (A3) 

where Uy is the number of years the experiment will run, A^^ is Avogadro's number and 
is the mass of the detector, in kilotons. 
Within this work we use: 

Uy = 5 years 
Mrf(z/^) = 50 Kt (A4) 

MdiUr) = 5 Kt 
2. Neutrino Flux, Oscillation Probability and Cross Section 

In order to simulate a realistic neutrino flux, we shall use the corresponding expressions 



for the z/g and z/^ unpolarized fluxes in the laboratory frame 
ddu,, An 
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Ely\l -/3 cos if) [3ml - 4^2^(1 - (3cosip)] (A5) 

^L^^fy^^ - /?cosy.) [ml - 2Ely{l - /^cosy.)] (A6) 

where is the number of useful muons, is the muon mass, L is the baseline, E^ is the 
parent muon energy, /3 = ^1 — mj^/ E'^, y = Ey/E^ and (p is the angle between the beam 
axis and the detector's direction. We assume that = 0° in the forward direction of the 
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muon beam, and consider an angular divergence 5(p = 0.1 mr. We also take = 2.5 x 10 
useful muons per year. The neutrino and antineutrino fluxes are identical. 

The neutrino oscillation probability is obtained by solving the Eq. ([5]) numerically, with- 



out any analytical approximations. We use the Preliminary Reference Earth Model 38| to 
simulate the electron density. 

The signal to be detected will take into account only charged current interactions. How- 
ever, the respective background for each signal can also involve neutral current interactions. 
In order to simulate all neutrino and antineutrino cross-sections, we interpolate numerical 



values for deep inelastic scattering found within 41|. 
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3. Neutrino Detection 



The choice of detector shall depend on the channel we are measuring. For the z/e and 
—>■ channels, we shall use a magnetized iron calorimeter, while for the z/g we will 
consider an emulsion cloud chamber. Nontheless, the simulation of both detectors is done 
in a similar way, by using Eq. (]A2|) . We take into account three main factors: the detection 
resolution Rp{Ey,E'^, the detection efficiency e^L^^) and the sources of background. 



Extensive studies have been made in 



30| and [l7| to describe z/^ and v.^ detection, respec- 



tively. A precise description of z/„ detection, including the background levels for v^. 



and z/,, 



z/^ channels, can be found in 



40[ |. which we follow. 



One defines the resolution function Rp{E^,El) in Eq. (]A2[) as a Gaussian distribution 
function: 



RfsiEiy, El^ 



exp 



'E\ 
2\l 



l\2 



(A7) 



where the effective relative energy resolution A/? depends on the detector, we assume it is 
the same for particles and antiparticles. For the z/g and z7^ channels, we take 

= O.lbE^, as in [40j. Regarding the channel resolution, we use the value found 

in [izl of = 0.20E^. 



The neutrino detection efficiency is taken from 401 ] and ITj. For the Ue — >■ and 
channels, we use the same functions as in ||40i] : 



K 

4 \ 4 



1 ) El<20 GeV 
E' > 20 GeV 



(A8) 



where we make a cut at energies lower than 4 GeV, and where is an efficiency normal- 
ization factor that corresponds to the high energy efficiencies presented in 40| for charged 
current interactions. Such factors appear in Table lllli 

n 

For the z/g —>■ Vj. efficiency, we use the data presented in |17l] to model a linear efficiency for 
charged current interactions. We scale our results in order to reproduce a global efficiency 
around 5%: 



6.69 X 10"^ {El + 20.86) 



(A9) 



where we conservatively assume that the same detection efficiency applies to neutrinos and 
antineutrinos. The behavior of the Vj. detection efficiency is shown in Figure [181 
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The background for each channel is calculated using Eq. (lA2p . replacing the efficiency 
with a background level. We present in Table [IV] the corresponding backgrounds for each 
channel, including their back grou nd level and a short description. More detail about these 
backgrounds can be found in SOj and 171 ]. 
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sin^ 2^23 


sin^ 26*13 


5 


• 


2a 


3<T 


fx^ • Y 


0.95 


0.1 


0° 
90° 


95.15 
92.29 


103.18 
100.32 


109.31 
106.45 


2200.73 
2900.53 


0.01 


0° 
90° 


96.17 
95.01 


104.20 
103.04 


110.33 
109.17 


226.83 
370.68 


0.99 


0.1 


0° 
90° 


92.91 
93.16 


104.20 
101.19 


110.33 
107.32 


2553.92 
3231.13 


0.01 


0° 
90° 


94.63 
97.84 


102.66 
105.87 


108.79 
112.00 


224.31 
389.24 



TABLE I: Minimum values obtained by assuming the correct (Xmin) ^^'^ wrong ((Xmm)"^) ^^S^ 
of Arri^i. The 2a and 3a confidence levels are defined with respect to the lowest x^. Notice that 
in the considered cases (Xmin)'" larger than the 3a level. 
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sin^ 2^23 


6 


reM 


= 1 


reM 


= 5 






2a 


3a 


2a 


3a 




0° 


0.015 


0.023 


0.005 


0.008 


0.95 


60° 


0.01 


0.025 


0.006 


0.01 




90° 


0.014 


0.033 


0.007 


0.01 




0° 


0.02 


0.05 


0.01 


0.015 


0.975 


60° 


0.02 


0.05 


0.01 


0.015 




90° 


0.025 


0.05 


0.01 


0.02 




0° 


0.06 


0.09 


0.02 


0.045 


0.99 


60° 


0.055 


0.095 


0.025 


0.07 




90° 


0.06 


0.105 


0.04 


0.07 



TABLE IL Minimum values of sin^ 2^13 with no degeneracies at a confidence level lower than 2a 
and 3(7, considering r^M = 1 and r^M = 5. 
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Channel 




I'e l^^i 


0.35 


Ve 


0.45 


l^iJi 


0.45 




0.35 



TABLE III: Normalization factors for charged current detection efficiency at a magnetic iron 
calorimeter {40]. 



Signal 


Background 




Description 


Background 


Channel 


Source 






Level 


l^e J^fi 


'^x 


NC 


Fake [i~ from ■k~ , K~ decay 


5 X 10-6 




Vf, 


cc 


No detection of /i+ and 
fake /x^ from hadron decay 


5 X 10-6 


Ve Vt 


Vfj, 


cc 


Induced charm production 


3.7 X 10-6 






cc 


Wrong charge identification 


1.8 X 10-'' 




I'x 


NC 


Hadron misidentification 


1 X 10-6 


v^i 


Vx 


NC 


Fake ^~ from vr", decay 


1 X 10"^ 



TABLE IV: Background levels considered for each signal channel 
levels and type of sources for antiparticles. 
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{e,,,S\ (0,„S), {9,„S), {0,„S), {d,„S), {9,„S), {9,„S), {0,„S) 

FIG. 1: Eight-fold degeneracy, as a product of three two-fold degeneracies. 




tan^023 

FIG. 2: Equiprobability curve for the Ve t'^ channel, showing the dependence of tan^ ^13 on 
tan^ ^23) for neutrinos (left) and antineutrinos (right). The probability considers sin^ 2^13 = 0.1 
(top) and 0.015 (bottom), which correspond to tan^ 613 ^ 0.026 and 3.8 x 10~^ respectively, with 
the other parameters at their 'true' values. We use L = 3000 km and Ei, = 30 GeV. 



28 



tan'0. 



FIG. 3: Same as Figure [2l but including the — > equiprobability curve (dashed). The in- 
tersection point indicates the real value of tan^ ^13 and tan^023) but notice that the value of b is 
required. 
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tan'0. 



FIG. 4: Intersection points of fe — > and fe Vt channels, for neutrinos (left) and antineutrinos 
(right), taking into account the variation of 5. The probability considers sin^ 2^13 = 0.1 (top) and 
0.015 (bottom), which correspond to tan^ ^13 ~ 0.026 and 3.8 x 10~^ respectively, with the other 
parameters at their 'true' values. We use L = 3000 km and E^, = 30 GeV. 
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FIG. 5: Intersection points for v^, — > and i/g v-r channels, for both signs of Am|-|^ at L = 
3000 km and Ey = 30 GeV. The probabihty considers sin^ 26*13 = 0.1 (left) and 0.015 (right), 
which correspond to tan^ ^13 ~ 0.026 and 3.8 x 10^'^ respectively, with the other parameters at 
their 'true' values. The red (green) line corresponds to neutrinos (antineutrinos) for the correct 
sign (+), while the blue (magenta) line corresponds to neutrinos (antineutrinos) for the wrong sign 
(— ) of Am|]^. We can see that the lines for the correct sign have the same value of tan^ ^13, while 
the lines for the wrong sign do not. 
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FIG. 6: Equiprobability curves for neutrinos (red) and antineutrinos (blue), for a fixed ^13, at 
L = 3000 km and Ei, = 20, 30, 40 GeV. We show equiprobability curves for the (green) and 

17^ (magenta), with Ey = 30 GeV. We have, from left to right, top to bottom, sin^ 2^13 = 0.1, 

0.05, 0.025 and 0.015, corresponding to tan^ 6*13 « 0.026, 0.013, 6.3 x 10"^ and 3.8 x lO^^. 
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FIG. 7: Degeneracies produced from the variation of ^^3. The probability is generated with 
sin^ 2^13 = 0.1, with the other parameters at their 'true' values, with L = 3000 km and 
Ey = 20, 30, 40 GeV. The values of sin^ 26'']^3 used to generate each set of curves are, from left 
to right, top to bottom, 0.1 (real value), 0.0668, 0.145 and 0.09 



31 



-100 100 -100 100 

6 

FIG. 8: Same as Figure [71 but with the probabihty generated by sin^ 2^13 = 0.05. The values of 
sin^ 26*^3 used are 0.05 (real value), 0.0341, 0.085 and 0.05. 
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FIG. 9: Same as Figure [TJ but with the probability generated by sin^ 2^13 = 0.025. The values of 
sin^ 26*^3 used are 0.025 (real value), 0.0176, 0.05 and 0.029. 
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FIG. 10: Same as Figure [3 but with the probabihty generated by sin^ 26*13 = 0.015. The values of 
sin^ 26*^3 used are 0.015 (real value), 0.01085, 0.035 and 0.0205. 




FIG. 11: Equiprobability curves for the 'magic baseline', in contrast to those for L = 3000 km 
(black). The curves in red (blue) correspond to neutrinos (antineutrinos) at L = 7250 km and 
= 20,30,40 GeV. We use, from left to right, top to bottom, sin^ 26113 = 0.1,0.05,0.025 and 
0.015. 
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FIG. 12: Expected allowed regions from the combination of data using three oscillation channels at 
2 (red) and 3a (blue). We use, from left to right, top to bottom, sin^ 26*13 = 0.01,0.05,0.025,0.015 
2.6 X 10-2,1.3 X 10-2,6.3 x 10-^,3.8 x 10-^). In this case we take r,M = 1- 
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FIG. 13: Same as Figure [T2\ but with r^Af = 3. 
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FIG. 14: Same as Figure [121 but with r^Af = 5. 
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FIG. 15: Sensitivity to solving the 623 degeneracy as a function of r^M- The red zones indicate 
regions where a degeneracy is present at a 2a level. The blue zones indicate regions where a 
degeneracy is present at a 3a level. We show results for sin^ 2^23 = 0.95 (top), 0.975 (middle) and 
0.99 (bottom), and 5 = 0° (left), 60° (center) and 90° (right). 
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FIG. 16: Same as Figure [T2t but including a Ue ^ i'^ measurement at the 'magic baseline', and 
using r,A/ = 1. 
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FIG. 18: Neutrino detection efficiency for Vt-. The average efficiency is under 5%, represented as 
eo- 
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